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ABSTRACT 

Context. We report on a detailed study of the Fe K emission/absorption complex in the nearby, bright Seyfert 1 galaxy Mrk 509. 
The study is part of an extensive XMM-Newton monitoring consisting of 10 pointings (~ 60 ks each) about once every four days, and 
includes also a reanalysis of previous XMM-Newton and Chandra observations. 
' Aims. We aim at understanding the origin and location of the Fe K emission and absorption regions. 

op I Methods. We combine the results of time-resolved spectral analysis on both short and long time- scales including model independent 

rms spectra. 

■ Results. Mrk 509 shows a clear (EW= 58 + 4 eV) neutral Fe Ka emission line that can be decomposed into a narrow (cr = 0.027 keV) 

' component (found in the Chandra HETG data) plus a resolved (cr = 0.22 keV) component. We find the first successful measurement of 

I a linear correlation between the intensity of the resolved line component and the 3-10 keV flux variations on time- scales of years down 

to a few days. The Fe Ka reverberates the hard X-ray continuum without any measurable lag, suggesting that the region producing 
the resolved Fe Ka component is located within a few light days-week (r < 10^ rg) from the Black Hole (BH). The lack of a redshifted 
wing in the line poses a lower limit of >40 rg for its distance from the BH. The Fe Ka could thus be emitted from the inner regions of 
_ the BLR, i.e. within the ~80 light days indicated by the 11/3 line measurements. In addition to these two neutral Fe Ka components, 

, we confirm the detection of weak (EW~ 8-20 eV) ionised Fe K emission. This ionised line can be modeled with either a blend of 

^ ■ two narrow Fe XXV and Fe XXVI emission lines (possibly produced by scattering from distant material) or with a single relativistic 

line produced, in an ionised disc, down to a few rg from the BH. In the latter interpretation, the presence of an ionised standard of-disc, 
down to a few rg, is consistent with the source high Eddington ratio. Finally, we observe a weakening/disappearing of the medium 
and high velocity high ionisation Fe K wind features found in previous XMM-Newton observations. 

Conclusions. This campaign has made possible the first reverberation measurement of the resolved component of the Fe Ka line, 
from which we can infer a location for the bulk of its emission at a distance of r~ 40 - 1000 rg from the BH. 

Key words, galaxies: individual: Mrk 509 - galaxies: active - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

X-ray observations of AGN have shown the almost ubiquitous 

presence of the Fe Ka line at 6.4 keV (Yaqoob et al. 2004; 

Send offprint requests to: Gabriele Ponti Nandra et al. 1997; 2007; Bianchi et al. 2007; de la Calle et al. 

e-mail: ponti@iasfbo.inaf. it 2010). Unlike the optical-UV lines that are emitted by distant 
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material only, the Fe Ka line traces reflection not only from dis- 
tant material (such as the inner wall of the molecular torus, the 
broad line region and/or the outer disc) but also from regions as 
close as a few rg (where rg=GM/c^) from the BH (Fabian et al. 
2000). 

The powerful reverberation mapping technique, routinely 
exploited on optical-UV lines (Clavel et al. 1991 ; Peterson 1993; 
Kaspi et al. 2000; Peterson et al. 2004), can also be applied to 
X-ray lines such as the Fe Ka line. This kind of analysis has a 
tremendous potential, allowing us to map the geometry of mat- 
ter surrounding the BH, starting from distances of a few gravita- 
tional radii up to light years. However, each Fe Ka component 
is expected to respond on a different characteristic time (years- 
decades for the torus, several days-months for the BLR-outer 
disc, and tens of seconds to a few hours for the inner accre- 
tion disc) and current X-ray instruments cannot easily disentan- 
gle the different components. Indeed, reverberation mapping of 
all Fe Ka emission components represents an enormous obser- 
vational challenge, and specially tailored monitoring campaigns 
(to sample the proper time scales) have to be designed. 

Since the detection of the first clear example of a broad and 
skewed Fe line profile in the spectrum of an AGN (indicating 
that most of the line emission is produced within few tens of rg; 
e.g. MCG-6-30-15, Tanaka et al. 1995) the quest to understand 
how the broad Fe Ka line varies with the continuum is ongo- 
ing. Indeed, close to the BH the simple one-to-one correlation 
between continuum and reflection line is distorted by General 
and Special relativistic effects. Several papers present extensive 
theoretical computations to describe the inner disc reverbera- 
tion to the continuum taking into account all relativistic effects 
(Reynolds et al. 1999; Fabian et al. 2000; Reynolds & Nowak 
2003). 

Several techniques have been employed to measure the 
variability-reverberation of the relativistic Fe Ka line. However, 
for the best cases such as MCG-6-30-15, the relativistic Fe line 
showed a complex behaviour, having a variable intensity at low 
fluxes (Ponti et al. 2004; Reynolds et al. 2004) while show- 
ing a constant intensity at higher fluxes (Vaughan et al. 2003; 
2004; see also the case of NG4051: Ponti et al. 2006). This 
puzzling and unexpected behaviour has been interpreted as due 
to strong light bending effects by some authors (Miniutti et al. 
2003; 2004) or, alternatively, as the evidence that the broad wing 
of the Fe Ka line is produced by strong and complex absorption 
effects (Miller et al. 2008). 

Thanks to the application of Fe Ka excess emission maps 
(Iwasawa et al. 2005; Dovciak et al. 2004; De Marco et al. 
2009), it has been possible to track weaker coherent patterns of 
Fe Ka variations. In a few sources Fe Ka variations are consis- 
tent with being produced by orbiting spots at a few r^ from the 
BH (Iwasawa et al. 2004; Turner et al. 2006; Petrucci et al. 2007; 
Tombesi et al. 2007). Future larger area telescopes are needed to 
finally assess if these features are present only sporadically dur- 
ing peculiar periods or if, instead, although weak, are always 
present and can be used to map the inner disc (see e.g. Vaughan 
et al. 2008; De Marco et al. 2009). 

A leap forward in X-ray reverberation studies occurred 
thanks to the application of pure timing techniques to the long 
XMM-Newton observation of 1H0707-495 that allowed the dis- 
covery of a "reverberation lag" between the direct X-ray contin- 
uum and the soft excess, probably dominated by FeL line emis- 
sion (Fabian et al. 2009; Zoghbi et al. 2010). Soon after simi- 
lar delays were seen in a few other objects (Ponti et al. 2010; 
De Marco 2011; Emmanouloulos et al. 2011; Zoghbi & Fabian 
2011; Turner et al. 2011). Recently, De Marco et al. (2012) 



showed that these lags are ubiquitous in AGN, that they scale 
with Mbh and have amplitudes of the order of the light crossing 
time of a few r^, thus suggesting a reverberation origin of the 
delay (but see also Miller et al. 2010). Another fundamental step 
forward will be to combine these timing techniques to detect re- 
verberation lags in the Fe K band (see Zoghbi et al. 2012). 

Reverberation from distant material has the advantage that 
the intensity of the Fe Ka line and the continuum are expected 
to follow a simple one-to-one correlation, however, the expected 
delays between the reflection component and the direct emis- 
sion are usually too large for a typical X-ray exposure. In fact, 
reflection from the inner walls of a molecular torus is expected 
to be delayed by a few years up to several decades and thus re- 
quires a very long monitoring campaign. Reflection from the 
BLR and/or outer disc is more accessible, the delay between 
continuum and reflection is expected to be between a few days 
up to few months. Thus a properly tailored monitoring campaign 
on a bright AGN with XMM-Newton, Chandra or Suzaku could 
achieve this goal. Several attempts have been made (Markowitz 
et al. 2003; Yaqoob et al. 2005; Liu et al. 2010). However, the 
15-20 % or larger error on the flux of the Fe Ka line and the low- 
sampling frequency of the X-ray observations have made the ap- 
plication of reverberation of the Fe Ka line on weeks-months 
timescales basically impossible, until now. 

Mrk 509 (z=0.034397) is one of the brightest Seyfert 1 
galaxies of the (2-100 keV) X-ray sky (Malizia et al. 1999; 
Revnivtsev et al. 2004; Sazonov et al. 2007), thus it has been ob- 
served by all major X-ray /Gamma-ray satellites. The Chandra 
HETG spectrum shows a narrow component of the Fe K line 
with an Equivalent Width (EW) of 50 eV (Yaqoob et al. 2004). 
XMM-Newton and Suzaku data provide evidence for a second 
broader (cr = 0.12 keV) neutral Fe K line (Ponti et al. 2009) as 
well as a weak ionized emission feature between 6.7-6.9 keV 
(Pounds et al. 2001; Page et al. 2003; Ponti et al. 2009). The 
ionised emission can be fit either using a relativistically broad- 
ened ionised line or an outflowing photo-ionised gas component. 

Imprinted on the Fe K band emission of Mrk 509 are the 
fingerprints of two kinds of ionised absorption components, 
one marginally consistent with a medium velocity outflow (v~ 
14000 km s"^ Ponti et al. 2009) and the others out(in)flowing 
with relativistic velocities (Cappi et al. 2009; Dadina et al. 2005; 
Tombesi et al. 2010). 

Here, we present the spectral and variability analysis of the 
Fe K complex energy band of Mrk 509 using the set of 10 XMM- 
Newton observations (60 ks each), about one every fours days 
and spanning more than 1 month, which we obtained in 2009 
(see the 3-10 keV light curve in Fig. [T]). We also re-analyse the 
previous 5 XMM-Newton observations. Thanks to this extensive 
monitoring campaign we can measure correlated variations be- 
tween the Fe K line intensity and X-ray continuum flux, allowing 
us, for the first time, to perform a reverberation mapping study 
on this X-ray emission line. In addition we can study the pres- 
ence of highly ionised matter from the innermost regions around 
the BH. 

The paper is organised as follows. §2 is devoted to the de- 
scription of the observations and data reduction. In §3 a first 
parametrisation (with a single Gaussian profile for the Fe Ka 
line) of the total summed spectrum of the 2009 campaign is pre- 
sented. Section 4 is dedicated to the detailed study of the Fe Ka 
emission. We first present the study of the Fe Ka line variabil- 
ity, assuming a single Gaussian profile (§4.1) we then use the 
Chandra HETG data (§4.2) to decompose the Fe Ka line in two 
Gaussian (narrow and resolved) components. §4.3 presents the 
correlation between Fe Ka intensity and the 3-10 keV contin- 
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Fig. 1. 3-10 keV EPIC-pn light curve (3 ks time bins) of the 10 
XMM-Newton observations of the 2009 monitoring campaign. 

uum (once the Fe Ka line is fitted with 2 Gaussian lines) which 
is confirmed, in a model independent way, by the rms spectrum 
(§4.4). In §4.5 we discuss the possible origin of the Fe Ka line. 
§5 presents the study and the discussion of the origin of the 
ionised Fe K emission/absorption. Conclusions are in §6. 

2. Observations and data reduction 

Mrk 509 was observed for a total of 15 times by XMM-Newton: 
on 2000-10-25, 2001-04-20, 2005-10-16, 2005-10-20, 2006- 
04-25, and 10 times in 2009 (see Fig. [H starting from 2009-10-15 
and ending on 2009-1 1-20). Ponti et al. (2009) and Kaastra et al. 
(2011) provide a full description of the first 5 and the last 10 
XMM-Newton observations, respectively. 

We initially reduced the EPIC data (as in Mehdipour et al. 
2011), starting from the ODE files, using the standard SAS v9.0 
software. However, we noted that the rest- frame best fit energy 
of the Fe Ka line in the EPIC-pn spectrum (EpeKa = 6.35 + 0.01 
keV) was not consistent with the best fit energy in the summed 
spectrum of the EPIC-MOS data (EpeKa = 6.41 ± 0.01 keV). 
This discrepancy (~ 50 eV) was found to be systematic, and was 
present in all 10 observations. Being significantly larger than the 
reported systematic uncertainty on the calibration of the absolute 
energy scale of 10 eV (CAL-TN-0018), this result triggered an 
in-depth study of the pn and MOS energy scales by the XMM- 
Newton EPIC calibration team. After excluding that this effect 
is related to X-ray loading, a stronger than expected long-term 
degradation/evolution of the Charge Transfer Inefficiency (CTI) 
was found. The pn long-term CTI was thus re-calibrated, and 
its corrected value implemented in the SAS version 10.0.0 (see 
CCF release note XMM-CCF-REL-271Q). 

The EPIC data were thus reduced again using the SAS 
version 10.0.0. During the XMM-Newton monitoring, both 
the EPIC-pn and the EPIC-MOS cameras were operating in 
the Small Window mode with the Thin filter applied. The 
#XMMEA_EP and #XMMEA_EM, for the pn and MOS cam- 
eras, respectively, are used to filter the events lists and to create 
Good Time Intervals (GTI). The FLAG==0 is then used for se- 
lection of events for making the spectra. The data were screened 
for increased flux of background particles. The contribution due 

^ S ee http : //xmm2 . esac . e sa. int/external/xmm_s w_cal/calib/rel_notes/ 
index, shtml 



to soft protons flares was negligible during the whole 2009 mon- 
itoring. The final cleaned EPIC-pn exposures for each XMM- 
Newton observation were about 60 ks, i.e. roughly 40 ks, after 
accounting for the proper dead-time of the pn when operating in 
small window mode (see Table 1 of Mehdipour et al. 201 1, for a 
list of the exposure times). 

The pn and MOS spectra were extracted from a circular re- 
gion of 45'' and 20'' radius centred on the source, respectively. 
The background was taken locally from identical circular re- 
gions located on the same CCD of the source for the EPIC- 
pn but on another CCD for the EPIC-MOS. The EPIC data 
showed no evidence of significant pile-up, thus single and dou- 
ble events were selected for both the pn (PATTERN<=4) and the 
MOS (PATTERN<=12) camera. Response matrices were gener- 
ated for each source spectrum using the SAS tasks arfgen and 
rmfgen. The sum of the spectra has been performed with the 
MATHPHA, ADDRMF and ADDARF tools within the HEASOFT 
package (version 6.10). 

Mrk509 was observed by the Chandra Advanced CCD 
Imaging Spectrometer (ACIS: Garmire et al. 2003) with the 
High-Energy Transmission Grating Spectrometer (HETGS: 
Canizares et al. 2005) in the focal plane, on 2001, April 13th 
(obsid 2087). Data were reduced with the Chandra Interactive 
Analysis of Observations (CIAO: Fruscione et al. 2006) 4.2 
and the Chandra Calibration DataBase (CALDB) 4.3.1 software, 
adopting standard procedures. 

All spectral fits were performed using the Xspec software 
(version 12.3.0) and include the neutral Galactic absorption 
(4.44x10^0 cm-2; Murphy et al. 1996), the energies are in the 
rest frame if not specified otherwise, however the energies in the 
plots are in the observed frame and the errors are reported at the 
90 per cent confidence level for one interesting parameter (Avni 
1976) in all the tables, while they are 1 a errors in the Figures. 
Mrk 509 has a cosmological redshift of 0.034397 (Huchra et al. 
1993) corresponding to a luminosity distance of 145 Mpc (tak- 
ing < ifo >= 73 km s"^ Mpc-\ Qa = 0.73 and = 0.27). 

3. The mean spectrum 

The upper panel of Figure [21 shows the data to best fit model ra- 
tio plot of the summed spectra of the 10 EPIC-pn observations 
performed during the 2009 XMM-Newton monitoring campaign, 
fitted in the 3.5-5 and 7.5-10 keV band with a simple power law, 
absorbed by Galactic material (interstellar neutral gas; phabs 
model in XSPEC). For comparison, the black data points in the 
lower panel of Fig. [2l show the same plot for the summed EPIC- 
pn spectrum of the previous 5 XMM-Newton observations taken 
between 2000 and 2006, while the red data show the summed 
XIS0+XIS3 spectra of the 4 Suzaku observations performed be- 
tween April and November 2006 (see Ponti et al. 2009 for more 
details). 

Thanks to a longer integrated exposure and a slightly higher 
flux, the source spectrum in the FeK band has significantly better 
statistics during the 2009 campaign than the sum of all the pre- 
vious observations (see Fig. O. Hence we can better constrain 
the FeK complex and study its variability not only on the time- 
scales of days and weeks over which the monitoring has been 
performed, but also on time-scales of years, using previous ob- 
servations. 

The upper panel of Fig. [2] shows, as observed during previous 
observations, an evident emission line at 6.4 keV as well as an 
emission tail at higher energies. A simple power law fit to the 4- 
10 keV total pn spectrum gives an un-acceptable fit (x^ = 2445.7 
for 1 197 dof). Adding a Gaussian line to fit the Fe Ka line at 6.4 
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Fig. 3. (Left, middle and right panels) Intensity vs. time in MJD, intensity vs. 3-10 keV flux and EW vs. 3-10 keV flux of the Fe Ka 
line fitted with a single Gaussian profile with cr = 0.092 keV (model 1), respectively. Each data point represents the best fit result 
obtained from the fit of the spectrum of each XMM-Newton observation (red hexagonals) and the total spectrum of the 2009 data 
(blue stars). The dashed lines represent the expected relations if the Fe Ka line has constant intensity, while the dash-dotted lines 
represent a constant EW, and the dotted line shows the best fit trend. The left panel show the observation date in Modified Julian 
Date, MJD minus 51500 which corresponds to 1999 November 18*. Two XMM-Newton observations occur at day 2159 and 2163 
and appear overlapping in the left panel. 



keV improves the fit by Ax^ = 1 134 for 3 extra parameters (see 
model 1 in Tab.[T]). This line reproduces the bulk of the neutral 
FeK emission (E = 6.43 ± 0.01 keV), but leaves strong residuals 
at higher energies which might be, at least in part, associated 
with an Fe KJ3 emission line. We thus add a second line with the 
same line width as the Fe Ka line and an energy of 7.06 keV 
(Kaastra & Mewe 1993). The best fit parameters are given as 
model 2 in Tab. [H and the fit improves by Ax^ = 66.9 for the 
addition of 1 new parameter. The best fit intensity is Npe = 
0.75 X 10"^ ph cm"^ s"\ and the observed Kj3/Ka ratio is = 
0.19 ± 0.02. This is slightly, but significantly, higher than 0.155- 
0.16, the value estimated by Molendi et al. (2003; who assume 
the Basko 1978 formulae; but see also Palmeri et al. 2003 a,b). 
The excess of Fe KJS emission might be due to the contamination 
by the Fe XXVI emission line at 6.966 keV. This idea is further 
supported by the observation of clear residuals, in the best fit, 
around 6.7 keV; and by the best fit FeK energy (E = 6.43 ± 
0.01 keV) which is inconsistent with the line arising from neutral 
iron. This thus suggests that the Gaussian line is trying to fit both 
the neutral and ionised Fe K components. 

To model ionised FeK emission, we add two narrow (cr = 0) 
emission lines, one (Fe XXV) emitting between 6.637 and 6.7 
keV (to take into account emission for each component of the 
triplet) and the other (Fe xxvi) emitting at 6.966 keV (model 
3 in Tab. [T]). Moreover we impose that the intensity of the KJ3 
has to be 0.155-0.16 times the intensity of Fe Ka one (and 
o'fq Kfi = crpe Ka)- The fit significantly improves (Ax^ = 25 for 
the addition of one more parameter). Both Fe XXV (the best 
fit line energy is consistent with each one of the triplet) and 
Fe XXVI are statistically required (although the Fe XXVI line 
is not resolved from the Fe KJ3 emission, thus its intensity de- 
pends on the assumed Ka/KjS ratio). In this model the Fe Ka line 
is roughly consistent with being produced by neutral or lowly 
ionised material (EpeKa = 6.415 ± 0.012 keV). 

Pounds et al. (2001), Ponti et al. (2009), de la Calle et al. 
(2010), Cerruti et al. (2011) and Noda et al. (2011) suggest that 
the inner accretion disc in Mrk509 might be highly ionised and 
thus the ionised emission of the Fe K complex might be asso- 
ciated to a relativistic ionised reflection component produced in 
the inner disc. To test this hypothesis, we substitute the two nar- 
row lines (Fe XXV and Fe xxvi) with one broad ionised line 
(model 4 in Tab. [TJ with a relativistic profile (disc line profile 



for a Schwarzschild black hole; diskline model in Xspec). We 
fixed the line energy either to 6.7 or 6.96 keV (for the Fe xxv or 
Fe XXVI line, respectively) and the inner and outer disc radius 
to 6 and 1000 gravitational radii. In both cases the best fit with 
this model suggests the inner accretion disc to be moderately in- 
clined ~ 33 - 18°, to have a fairly standard disc emissivity index 
-2.2 - 2.8 and an equivalent width of the line to be EW~ 36 - 40 
eV. The data are described reasonably well by this model result- 
ing insix^ = 1216.0 and;^^ = 1219.4 for 1191 dof, for the 
case of a broad Fe xxv and Fe xxvi line, respectively. Thus, 
the single broad ionised disc-line model is statistically indistin- 
guishable from the multiple narrow lines one (x^ = 1219.8 for 
1192 dof). 

3.1. Medium velocity-high ionisation winds 

The summed spectrum of the previous XMM-Newton obser- 
vations showed a medium velocity (Vout ~ 0.048 ± 0.013 c) 
highly ionised (Log(^)~ 5) outflow in Mrk 509 (Ponti et al. 
2009). The associated Fe xxvi absorption line was detected 
both in the EPIC-pn and MOS camera, with an equivalent width 
EW= -13.1^2 9 with a total significance between 3 - 4(T 

(~ 99.9 % probability). 

During the 2009 XMM-Newton campaign this highly ionised 
absorption component is not significantly detected. If we add a 
narrow Gaussian absorption line at 7.3 keV, the energy of the ab- 
sorption feature in the previous XMM-Newton observations, we 
observe the line to be much weaker with the best fit line EW be- 
ing -3.2^2 g eV, significantly smaller than observed in previous 
XMM-Newton observations. 



4. The neutral Fe Ka component 

In this section we further investigate the nature of the Fe Ka line, 
looking at the individual spectra obtained over the years. Our 
analysis of the long 2009 monitoring campaign, which triples the 
total exposure on Mrk509, confirms the presence of a resolved 
component (cr = 0.092 ±0.012 keV) of the neutral Fe Ka line 
(see model 3, but also 4 and 5 of Tab.[T]). 
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Fig. 4. (Left and right panel) The intensity and EW of the resolved (cr = 0.22 keV) Fe Ka component (the narrow Fe Ka component 
is assumed to have cr = 0.027 keV and be constant) as a function of the 3-10 keV flux. Blue stars and red hexagonals show the best 
fit results of the 10 observations of the 2009 campaign and of the previous XMM-Newton observations, respectively. Dashed and dot 
dashed lines show the constant intensity and constant EW cases, respectively. The best fit relation (dotted line) is consistent with 
the resolved component of the Fe Ka line having constant EW. The intensity of the resolved component of the Fe Ka line follows 
the 3-10 keV continuum variations with a 1-to-l relation. 



4.1. Fe Ka variations on years time-scales 

To study the variability of the neutral Fe Ka emission line (on 
years time- scales), we fit the spectrum of each of the old XMM- 
Newton observations (i.e. between 2000 and 2006) with a single 
FeKof (plus associated KJ3 emission) plus two narrow emission 
lines (such as in ^in model 3) to parametrize the ionised FeK 
emission (see Tab.[T]). Leaving the Fe Ka widths free to vary, as 
in Tab. [H would result in unconstrained values for the spectra 
with the shortest exposures (due to the lower statistics). We thus 
decide to fix the width of the Fe Ka line to cr = 0.092 keV, 
its best fit value as observed in the total spectrum of the 2009 
campaign (model 3 of Tab. [T]). 

The left panel of Fig. [3] shows the Fe Ka line intensity for 
each XMM-Newton observation as a function of time, using the 
summed 2009 data. The line intensity is observed to vary by less 
than 25 %. The fit with a constant Fe Ka intensity (dashed line), 
however, is unsatisfactory (x^ = 10.7 for 5 dof). The middle 
panel of Fig. [3] shows the Fe Ka intensity vs. source flux in the 
3-10 keV band. The fit slightly improves when a linear relation 
(see dotted line) is considered (Ax^ = 7.6 for the addition of 
one more parameter; 97 % F-test probability). The increase in 
Fe Ka intensity with flux might suggest that the line is respond- 
ing quickly to the illuminating continuum, keeping a constant 
EW with flux. The dash-dotted line shows the expected Fe Ka 
intensity variation for a line with constant EW. The observed line 
intensity variations are intermediate between the constant inten- 
sity and constant EW cases. The right panel of Fig. [3] confirms 
that the line has neither a constant intensity nor constant EW, 
instead it sits somewhere in the middle between these two cases. 

The Fe Ka variations on years time- scales suggest that at 
least part of the line is varying following the 3-10 keV contin- 
uum. We want to point out that the width of the Fe Ka line is 
comparable to the EPIC-pn energy resolution. This means that 
the observed Fe Ka variability may be the product of a con- 
stant narrow component, coming from distant material, plus a 
broader, resolved and variable Fe Ka line produced closer to the 
BH. Unfortunately, due to the limited energy resolution of the 
EPIC cameras aboard XMM-A^^w^(9^, we cannot resolve the Fe K 
emission, coming from regions located at light weeks from those 



at light years from the BH, based on the Fe Ka line widths. Only 
with the Chandra high energy transmission grating (HETG) we 
can confidently pose some constraints on the distance of the dif- 
ferent Fe K emission components. 



4.2. Chandra HETG 

Chandra observed Mrk 509 with the HETG instrument only 
once for 50 ks (Yaqoob et al. 2003). During the HETG obser- 
vation the 3-10 keV flux was 4.36 x 10"^^ erg cm"^ s"^ with 
a power law spectrum of index F = 1.76^q q2. An excess was 
present at 6.4 keV, thus we added a Gaussian lin^l. In agreement 
with the results obtained by Shu, Yaqoob & Wang (2010) and 
Yaqoob & Padmanabhan (2004) we detect a line at 6.42 ± 0.02 
keV with an intensity of 3 ± 2 x 10"^ photons cm"^ s"^ . The line 
is resolved and has a width significantly smaller than the one 
measured by XMM-Newton, cr = 0.027+^-^|^ keV. This suggests 
that, at least part of, the neutral Fe Ka emission is produced in 
regions more distant than a few thousands gravitational radii. In 
fact, if the material is in Keplerian motion and assuming a BH 
mass of Mrk 509 of Mbh = 1.4 - 3 x 10^ M© (Peterson et al. 
2004; Mehdipour et al. 2011), then the narrow core of the line 
is produced at a distance of r=0.2-0.5 pc (~ 30000 r^). We note 
that this value is of the same order of magnitude as the molec- 
ular sublimation radius for Mrk 509. Landt et al. (2011) using 
quasi- simultaneous near-infrared and optical spectroscopy esti- 
mate a radius of the hot dust of ~ 0.27 pc (0.84 ly), which is also 
consistent with the one estimated following eq. 5 of Barvainis 
(1987) assuming a bolometric luminosity Lboi = 1.07 x 10"^^ erg 
s"^ (Woo & Urry 2002). This suggests that this narrow compo- 
nent of the Fe K line might be associated to the inner wall of the 
molecular torus. 



^ In order not to lose the excellent energy resolution, we decided not 
to rebin the spectrum and to fit the data using the C- statistics (Cash 
1979). The best fit has C-stat value of 370.5 for 284 dof. 
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Fig. 2. (Upper panel) Summed EPIC-pn spectra of the 10 ob- 
servations performed during the 2009 XMM-Newton monitoring 
campaign. The data are fit, in the 4-5 and 7.5-10 keV bands, with 
a simple power law, absorbed by Galactic material, and the ratio 
of the data to the best-fit model is shown. The dotted lines indi- 
cate the rest frame positions of the Fe Ka, Fe KJS, Fe XXV (reso- 
nance, intercombination and forbidden) and Fe XXVI lines. The 
X-axis reports the observed-frame energy. (Lower panel) In red 
the 2006 XIS0+XIS3 Suzaku summed mean spectra are shown. 
In black the summed spectra of the 5 XMM-Newton observa- 
tions performed between 2000 and 2006 are shown. This Figure 
is taken from Ponti et al. (2009). The arrows mark possible ab- 
sorption features in the spectrum. 



4.3. Two components of the Fe Ka line 

Thus, as suggested by the analysis of the Chandra HETG data 
and in agreement with the observed variability on years time- 
scales, we interpret the Fe Kof line as being composed by two 
components which are indistinguishable at the EPIC resolution. 
We first re-fit the mean spectrum of the 2009 campaign with two 
components for the Fe Ka line. One "narrow" component with 
the line width fixed at the best fit value derived from Chandra 
analysis (cr = 0.027 keV) plus a "resolved" component with 
its width free to vary. Model 6 in Tab. [21 shows the best fit re- 
sults assuming the energy of both components is the same. The 
narrow component has a best fit EW= 27 + 4 eV, while the re- 
solved neutral line has an EW= 42^^ eV and a best fit line width 



(T = 0.22 ± 0.05 keV which is larger than the single Gaussian 
Fe Ka fit (cr = 0.092 keV). 

Reflection is the most probable origin of the Fe Ka line. 
Associated to reflection lines an underlying reflection contin- 
uum is expected and generally observed. In particular the ra- 
tio between the intensity of the line over the reflection contin- 
uum strongly depends on the reflector column density reaching a 
value of EWpeKor ~ 1 keV for Compton thick materials. To check 
the impact of the reflection continuum on the best fit model we 
add a standard neutral reflection continuum (PEXRAV in XSPEC) 
with intensity such that the EWpeKo^ = 1 keV over their reflec- 
tion continua. The new best fit line EW do not vary significantly 
(EW= 26 ± 3 eV and EW= 33^^ eV, for the narrow and resolved 
Fe Ka lines, respectively). Thus, and considering also the lim- 
ited energy band used here we decide to disregard the continuum 
reflection component (the impact of the reflection component on 
the broad band source emission will be studied by Petrucci et al. 
2012 taking into account the UV to soft Gamma ray emission 
with physical models). The narrow component, as observed dur- 
ing the 2009 campaign, has an intensity of 1.5 ± 0.2 x 10"^ ph 
cm"^ s"\ which is about half the total Fe Ka intensity (see Tab. 
[T]and Fig.O. This value is consistent with the intensity of the 
narrow Fe Ka line observed by Chandra. The width of the nar- 
row component suggests a distance of 0.2-0.5 pc from the BH, 
thus we expect that all the variability on time- scales shorter than 
a few years would be smeared out because of light travelling ef- 
fects. For this reason and because the lower signal to noise in 
the individual 2009 and earlier XMM-Newton spectra does not 
allow us to disentangle both components, in all the following fits 
we assume a constant intensity of 1.5 x 10"^ cm"^ s"^ for this 
narrow Fe Ka component. 

Next, we fit the spectrum of each of the 15 XMM-Newton 
observations with model 6 shown in Tab. [21 assuming a constant 
intensity for the narrow Fe Ka component and constant width 
for the resolved component, plus the associated Fe Kfi and nar- 
row Fe XXV and Fe xxvi lines. The left panel of Fig. [H shows 
the intensity of the resolved component (cr = 0.22 keV) vs. the 
3-10 keV flux for the 10 observations of the 2009 monitoring 
(blue stars) as well as for the previous XMM-Newton observa- 
tions (red hexagonals). The dashed line shows the best fit as- 
suming that the Fe Ka intensity is constant, which results in an 
unsatisfactory fit with;^^ = 33.2 for 14 dof. On the other hand, 
once the data are fitted with a linear relation (dotted line), the fit 
significantly improves (A;^^ = 23.0 for the addition of 1 new pa- 
rameter, which corresponds to F-test probability > 99.98 %). We 
also compute the Pearson's linear correlation coefficient is equal 
to 0.87 and has a probability = 2.7 x 10"^, which corresponds to 
a significance of the correlation of more than 4 cr (similar results 
are obtained using a Spearman's rho or Kendall's tau correlation 
coefficients). The slope of the observed best fit relation is consis- 
tent with that expected if the resolved line is responding linearly 
to the continuum variations. This is confirmed in the right panel 
of Fig. [4l which shows the resolved Fe Ka line EW is consis- 
tent with being constant (x^ = 11. 1 for 14 dof), as expected if 
the Fe Ka line is responding linearly to the 3-10 keV continuum 
flux variations. 

The line intensity is significantly variable even on time- 
scales of few days, e.g. between the different pointing of the 
2009 monitoring campaign. In fact, fitting the 2009 Fe Ka in- 
tensities with a constant gives slx^ = 13.7 for 9 dof, which be- 
comes = 5.3 when a linear relation is considered (as specified 
in §2, conservative 90 % errors are used here). The Pearson's lin- 
ear correlation coefficient results to be 0.8 and has a probability 
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4 6 
Observation number 

Fig. 5. Intensity vs. observation number of the 10 XMM-Newton 
pointings of the 2009 campaign. The black dashed line shows 
the re-scaled (with the mean 3-10 keV flux equalling the mean 
line intensity) 3-10 keV source flux. The intensity of the resolved 
component of the Fe Ka line follows, with a 1-to-l relation the 
3-10 keV continuum variations without any measurable lag. Icr 
errors are shown. 



= 5 X 10"^, which corresponds to a significance of the correlation 
of about 3 cr. 

Fig. [5] shows the variations of the intensity of the resolved 
component of the Fe Ka line as a function of time (I cr errors 
are shown here), during the 2009 campaign, overplotted on the 3- 
10 keV rescaled flux (dashed line). As already suggested above 
and Fig. (H the Fe K line variations track that of the continuum 
very well. We also note that no measurable lag is present, thus 
this broad component of the Fe Ka line responds to the X-ray 
continuum within less than four days. 

4.4. Total RMS spectra 

Figure [6] shows the total Root Mean Square variability (rms) 
spectrum calculated between the 10 different observations of the 
2009 campaign. The rms has been calculated with ten time bins, 
each one being a 60 ks XMM-Newton pointing. Thus this rms is 
sampling the variability within the observation separation time- 
scale of about 4 days and the monitoring time-scale of slightly 
more than one month (see Fig.[T]). The total rms shows the spec- 
trum of the variable component, only. Thus, in contrast to the 
mean spectrum, it has no contribution from the constant emis- 
sion from distant material (i.e. the narrow core of the Fe K line). 
The uncertainties on the total rms are derived from the uncer- 
tainties on the fractional variability (see formula B. 2 of Vaughan 
et al. 2003; A.l of Ponti et al. 2004) multiplying for the mean 
and taking into account its error. 

The 3-10 keV total rms spectrum has a power law shape 
with spectral index F = 1.98 ± 0.06 and normalisation of 
1.05 ± 0.05 X 10"^ ph. cm"^ s"^ (see Fig. [6]). A clear excess 
of variability is present at 6.4 keV. The addition of a Gaussian 
line significantly improves the fit (Ax^ = 10.6 for the addition of 
2 dof, that corresponds to a F-test significance of 99.3 %). The 
best fit energy of the line is £" = 6.45 ± 0.08 keV and intensity 
(2 ± 1) X 10"^ ph. cm"^ s"^ The width of the line is constrained 
to be less than ~ 0.2 keV, consistent with the variable resolved 
neutral Fe K line. The equivalent width in the total rms spec- 
trum is EW = 71 ± 36 eV consistent with the EW of the resolved 
component observed in the mean spectrum. The detection of this 
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Fig. 6. Long time- scales total rms spectrum between the differ- 
ent observations. The spectrum is calculated with time bins of 60 
ks (corresponding to the exposure of each XMM-Newton point- 
ing) and the total monitoring time of about 1 month. An excess 
of variability is clearly evident at 6.4 keV (EW= 71 + 36 eV). 
This excess confirms, in a model independent way, the correlated 
variability of the resolved Fe K line component (see Fig.O. 



excess of variability indicates, in a model independent way, that 
the resolved component of the line is varying linearly with the 
continuum on these time-scales. This reinforces the robustness 
of the detection of Fe K reverberation. 



4.5. Locating the Fe Ka emitting region 

Are the observed variability properties in agreement with the 
spectral ones? Assuming that the material producing the re- 
solved Fe K line is in Keplerian motion around the BH, the line 
width (cr = 0.22 keV) impHes that it is located at 300 - 1000 
r^ from the BH. Assuming a BH mass of Mrk 509 of Mbh = 
1.4 - 3 X 10^ Mo (Peterson et al. 2004; Mehdipour et al. 2011), 
then this distance corresponds to about few days - a light week. 
The typical spacing between the different XMM-Newton obser- 
vations during the 2009 monitoring is about 4 days, thus it is 
in very good agreement with the observed Fe Ka variations. 
Moreover, the fast response of the Fe Ka flux to the continuum 
changes indicates that the bulk of the resolved Fe Ka emission is 
produced around or within several hundreds up to few thousands 
gravitational radii from the central BH. 

It is more difficult, however, to pose a lower limit to the po- 
sition of the Fe Ka emitting region. We note that if the Fe Ka 
line emitting region extends down to few gravitational radii from 
the BH, then the line shape should present a prominent red wing. 
Thus, we fit the resolved component of the Fe Ka line with a disc 
line profile (DIS KLINE in Xspec). In the fitting process we allow 
the line energy to vary in the range 6.4-6.42 keV, we fix the disc 
outer radius to 1000 r^ and the illumination profile to a = -3 (the 
expected value for a standard alpha disc; Laor 1991; Wilkins & 
Fabian 201 1). The lack of relativistic redshifted Fe Ka emission 
suggests an inner radius larger than ~ 85, 45 and 37 rg (which 
corresponds to roughly 8-16 light hours) for a disc inclination of 
30, 20 and 10 degrees, respectively. Thus, the resolved compo- 
nent of the Fe Ka line is probably emitted between 40 - 1000 r^ 
from the BH. 
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Can such a "narrow" disc annulus produce an Fe Ka line 
of ~ 40 - 50 eV? Reflection from an accretion disc with Solar 
iron abundances covering half of the sky is expected to pro- 
duce a Fe Ka line with EW~ 100 - 150 eV (Matt et al. 
1991). The line EW is expected to decrease/increase roughly 
linearly/logarithmically for iron abundances lower/higher than 
Solar (Matt et al. 1996; 1997). Steenbrugge et al. (2011) mea- 
sured in Mrk 509 a relative iron to oxygen abundance of Fe/0= 
0.85 ± 0.06. Assuming that this translate into an iron abundance 
of 0.85 Solar (but see Arav et al. 2007), this would correspond 
to EWpeKa ~ 90 - 130 eV. If the primary X-ray source in Mrk 
509 is compact (as the variability suggests; McHardy et al. 2006; 
Ponti et al. 2012) and located at a few r^ above the BH and if the 
disc is flat, we can estimate (neglecting relativistic effects) the 
geometric solid angle covered by the disc annulus producing the 
Fe Ka line (rin ~ 40 and rout ~ 1000 r^). If the primary X-ray 
source is located between 1 and 4 r^ above the BH (De Marco et 
al. 2012), the flat disc annulus covering factor would be between 
2 and 5 % of the sky. Thus reflection from such a flat annulus 
would produce (even in the extreme case of a Fe Ka EW= 130 
eV for a standard disc) a line with EW~ 4 - 13 eV. The observed 
EW of the resolved- variable Fe Ka line is several times larger 
(EW=42 eV) than this estimated value. This indicates a larger 
covering factor of the reflector, compared to the flat disc, sug- 
gesting that the material producing the Fe Ka line is distributed 
azimuthally above the disc, possibly in the form of clouds, per- 
haps associated to the inner BLR (see Costantini et al. 2012 for 
more details). 

The observed correlation on days-weeks timescales con- 
straints also in which part of the BLR the Fe Ka line is produced. 
We can exclude, in fact, that the Fe Ka emission is produced in 
the optical BLR (producing the bulk of H/S emission), because 
the HyS line is observed to reverberate with a delay of 80 days, 
being thus significantly more distant than the region producing 
the Fe Ka line. On the other hand, several studies show that the 
BLR might be stratified, with the higher ionisation lines located 
closer to the central BH. The correspondence between Fe Ka 
and the inner BLR is reinforced by the consistency between the 
width of the Fe Ka line (cr = 0.21 ±0.07 keV, which corresponds 
to FWHM ~ 1.5-3 X 10"^ km s"^) and the ones of the broadest 
components of the UV broad emission lines (e.g. hya, C IV, C 
III and O VI) which have components with FWHM ~ 10^ km 
s-i (Kriss et al. 2011). 

The lack of relativistic effects on the shape of the Fe Ka line 
suggests the absence of a neutral standard thin accretion disc ex- 
tending down to a few gravitational radii from the BH. However, 
this appears to be at odds with the high efficiency (Lboi -5-10 
% LEdd) of the disc emission of Mrk 509 (Mehdipour et al. 201 1 ; 
Petrucci et al. 2012). This leads to the question of why we do not 
find traces of the inner accretion disc in the Fe Ka line shape if 
it is present in this source. 

5. The ionised Fe K emission 

During both the 2009 campaign and the previous XMM-Newton 
observations, Mrk 509 clearly showed an excess of emission 
around 6.7-7 keV (see Fig. [2]) most probably associated to emis- 
sion from ionised iron. As shown in ^this excess can be mod- 
elled both by the combination of narrow emission lines from 
Fe XXV and Fe XXVI, or by a single relativistic emission line 
(see Tab. [TJ. The parameters of this weak ionised emission 
line(s) can be affected by the modelling of the stronger Fe Ka 
line. For this reason, we now re-fit the mean spectrum including 
both the narrow and the resolved component of the Fe Ka line. 



We first consider that the ionised emission is produced by 
narrow emission lines (Fe XXV and Fe XX Vl). Such emis- 
sion lines from highly ionised ions are now often observed 
(Costantini et al. 2010; e.g. for a compilation of sources, Bianchi 
et al. 2009a,b; Fukazawa et al. 2011) and they can arise from 
photo-ionised (Bianchi et al. 2005; Bianchi & Matt 2002) or 
collisionally ionised plasma (Cappi et al. 1999). Thus we fit the 
spectrum with two components for the Fe Ka line and the as- 
sociated Kfi lines, plus two narrow (cr = I eV) Gaussian emis- 
sion lines one (Fe XXV) with energy constrained to be between 
6.637 and 6.7 keV and the other (Fe xxvi) with energy fixed at 
E = 6.966 keV (see model 6 in Tab. [2]). The model reproduces 
the data wefl (x^ = 1198.5 for 1190 dof). The weakness of the 
Fe XXV and Fe XXVI lines prevents us from significantly con- 
straining the line variability between the different XMM-Newton 
observations. 

For comparison, we fit with the same model also to the 
summed spectrum of all XMM-Newton observations taken be- 
tween 2000 and 2006 (see model 7 of Tab.O. The ionised emis- 
sion lines are consistent with being constant within the two sets 
of observations. However, the statistics is not good enough to 
discriminate if it is the line intensity (which would suggest an 
origin at large distances) or the EWs remain constant. Clear 
is, instead, the variation of the medium outflow velocity highly 
ionised absorption line, which almost disappeared during the 
2009 campaign. The addition of this component to the model 
used to fit the combined spectra of the 2009 campaign, improves 
the fit by Ax^ =4.5 for 2 new parameters which corresponds to 
a F-test probability of ~ 90 %. Thus we decide to disregard this 
absorption component in all subsequent fits. 

Another clear difference compared to previous observations 
is related to the disappearance of the highly ionised absorp- 
tion with mildly relativistic (up to 0.14-0. 2c) outflow velocities 
(Dadina et al. 2005; Cappi et al. 2009; Tombesi et al. 2010). 
We searched, in fact, for such features in all 10 observations 
obtained after the XMM-Newton campaign, by including nar- 
row absorption lines in the model between 4-10 keV. We found 
only a marginal (Ax'^ ~6) detection of two absorption features 
at 9 keV and 10.2 keV (rest- frame energies) during observation 
4. Even if consistent with being produced by Fe XXVI Ka and 
KJS at v~0.3c, and similar to earlier results (Cappi et al. 2009), 
the level of (highly) ionised absorption during the XMM-Newton 
campaign is found to be significantly reduced compared with 
most previous XMM-Newton observations. We obtained upper 
limits (at 90 % confidence) on the equivalent width of narrow (cr 
fixed to 100 eV) Gaussian absorption lines with typical values 
between -5 and -30 eV, between 7.5 and 9.0 keV, depending on 
the energy and observation considered. This is typically lower 
than values (between -20 and -30 eV) found in the lines detected 
in earlier observations (Cappi et al. 2009, Tombesi et al. 2010) 
excluding that such UFOs were present during the 2009 cam- 
paign. 

5. 1. Collisionally ionised plasma 

We attempt to interpret the highly ionised emission lines via a 
self-consistent physical model. First we applied the collision- 
ally ionised model CI^ in spex (Kaastra et al. 1996). In this 
fit we considered the 3.5-10 keV band for the continuum. We 
used Gaussian components for the Fe Ka line profile and Fe 
K/3, constraining the flux of the latter to be 0.155-0.16 times 
the Fe Ka one (Palmeri et al. 2003). This was done in order 

^ http://www. sron.nl/files/HEA/SPEX/manuals/manual.pdf 
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to mitigate the degeneracy induced by the partial blend with the 
Fe XXVI Ly a line. The best fit points to a high-temperature gas 
(kT= 8.5 + 1 .5 keV). At this temperature, the predicted line fluxes 
are ~ 5 X 10"^phcm"^ s"^ and ~ 3 x 10"^phcm"^ s"^ for the 
Fe XXV triplet and the Fe XXVI Ly a, respectively. These values 
are consistent with those measured empirically using Gaussian 
lines (i.e. Table O. In theory these lines may be produced by hot, 
line emitting, gas in the form of a starburst driven wind. Mrk 509 
has a total luminosity L2-10 keV = 1-3 x 10"^"^ erg s"^ in the 2-10 
keV band. Assuming an Fe abundance of 0.4 Solar (as observed 
in starburst galaxies, Cappi et al. 1999) the best fit thermal star- 
burst model requires a luminosity L2-10 keV = 3.3 x 10"^^ erg s"^ 
to reproduce the Fe xxv and Fe XXVI line emission (reducing to 
L2-10 keV = 1.6 X 10"^^ erg s"^ for Solar iron abundance). Using 
the correlation between L2-10 keV and the far infrared luminosity 
(Lfir), valid in star forming galaxies (Ranalli et al. 2003), we 
estimate a corresponding Lfir > 10"^^ erg s"^ and a star forma- 
tion rate higher than 400 M© yr"\ which is several times larger 
than the actual total IR luminosity of Mrk 509, Ljr ~ 2 x 10^^ 
Lo ~ 8 X 10"^"^ erg s"^ (Rieke et al. 1978), thus we disfavour this 
interpretation. 

5.2. Photo-ionised plasma 

Alternatively, the highly ionised lines may be produced by a pho- 
toionised plasma. To test this we used a grid of parameters cre- 
ated using Cloudy (iFerland et al.lll998]) where the column den- 
sity logA^H of the gas ranged between 21.7 - 24.5 cm"^ and the 
ionisation parameter log(^) ranged between 3.4-7. The grid has 
been calculated using a covering factor of one. Since the intrinsic 
line luminosity scales linearly with the covering factor, we used 
the ratio between the model and the data of the Fe xxvi line as a 
reference for the covering factor. In Mrk 509, only Fe XXVI and 
Fe XXV have significant detection, while for other narrow lines 
from highly ionised ions (e.g. Ovill at 18.97 A and NeX at 
12.13 A) we have obtained upper limits from the RGS spectrum. 
These limits are useful in constraining the model (i.e. Costantini 
et al. 2010). In Fig. [7] we compare the line luminosities observed 
with those computed for a range of models which can fit the data. 
In order to reproduce the luminosity of the highly ionised iron 
ions, the gas should have log(^)= 4-5.1 and A^h = 23.4 - 24.2. 
The covering factor is Cy = 0.3 -0.5. As we do not see any asso- 
ciated absorption, the gas must be out of the line of sight. Such 
lines might possibly originate in e.g. the narrow line region or 
the highly ionised skin of the torus (Bianchi et al. 2005; Bianchi 
& Matt 2002). 



5.3. Ionised reflection from inner disc 

As shown in ^ the ionised emission can be fitted equally well 
with an ionised relativistic emission line. A comparably good fit 
is also obtained with a broad Gaussian Fe XXV profile in addi- 
tion to the double Fe Ka-\-p lines. We fix the energy of the broad 
Fe XXV line to E= 6.7 keV. The line is significantly broadened 
cr = 0.23 keV and moderately intense (EW= 15 eV). Although 
this broadening is not as extreme as to exclude a simple Compton 
broadening on the ionised surface of the accretion disc, we de- 
cided, as an alternative to the Gaussian line, to fit the ionised 
emission with a relativistic disc line profile. An acceptable fit is 
obtained also in this case (x^ = 1204.7 for 1 192 dof). Assuming 
a standard disc emissivity index, inclination and outer disc ra- 
dius ofjS = -3, a = 30° and rout = 200 r^, the best fit disc inner 




Fig. 7. Photoionization modeling for the highly ionized iron ions 
Fexxvi and Fexxv. Triangles: data. Light shaded line: range 
of models viable to fit the data. 

radius is rin = 21 Vg, consistent with a value as small as rin = 7 r^ 
and the disc-line equivalent width is, EW = 20^^ eV. 

It is difficult to estimate the line EW expected from an 
ionised inner annulus of the disc. In fact the Fe XXV and Fe XXVI 
line EW strongly depend on the poorly constrained disc ion- 
isation parameter (Garcia et al. 2011) and on the disc annu- 
lus covering angle. However for reasonable values of these pa- 
rameters (annulus covering angle ~ 20 - 40 % of the sky and 
log(^)~ 2.8 - 3.5 erg cm"^ s"\ which corresponds to the peak 
of Fe XXV and Fe xxvi emission) the line EW is expected to be 
between 5-50 eV. These results suggest that the inner part of the 
accretion disc might be highly ionised, thus explaining the lack 
of detection of relativistic Fe Ka line, combined with the high 
source efficiency (which suggests a thin standard accretion disc 
extending all the way down to the last stable orbit). 

6. Discussion and Conclusions 

We investigated the spectral variability of the Fe K band in the 
nearby, bright Seyfert 1 galaxy Mrk 509, using the 10 observa- 
tions of the 2009 XMM-Newton monitoring campaign as well as 
all the previous XMM-Newton observations, totalling an expo- 
sure of more than 900 ks in about 10 years, resulting in one of 
the best quality Fe K spectra ever taken of a Seyfert 1 galaxy. 
This allows us, for the first time, to perform reverberation map- 
ping of the resolved Fe Ka line. 

Figure [8] summarise in a sketch a possible scenario for the 
production of the Fe K emission in Mrk 509. The width of the 
narrow core of the Fe Ka line suggests an origin from distant 
material, possibly the inner wall of the molecular torus located 
at 0.2-few pc. The correlated variations (on a few days time- 
scales) between the 3-10 keV continuum and the intensity of the 
resolved component of the Fe Ka suggest an origin between sev- 
eral tens and few thousands of rg from the BH. The resolved Fe 
Ka emission can be produced in the disc, but we favour an origin 
at the base of a stratified broad line region. We note that none of 
the X-ray or UV absorption components with measured location 
is CO- spatial with the resolved Fe Ka emitting region. Moreover, 
the properties of the X-ray and UV absorbers appear to differ 
from the ones required to produce the resolved Fe Ka line, sug- 
gesting that this emitting material is outside the line of sight. 
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Fig. 8. Sketch of a possible locations for the different regions 
producing Fe K emission (diagram not to scale). The star repre- 
sents the primary X-ray source, located close to the BH. 

possibly in the form of an equatorial disc flattened wind such as 
observed in stellar mass black holes in the soft state (Ponti et. al. 
2012) and neutron stars (Diaz-Trigo et al. 2006). The ionised Fe 
K emission might be produced either by photo-ionisation from 
distant material, such as the narrow line region and/or the ionised 
skin of the torus, or in the ionised inner accretion disc. 
The results of this study show that: 

- The XMM-Newton spectrum of Mrk 509 shows an evident 
Fe Ka line with total EW= 58 ± 4 eV. Fitted with a single 
Gaussian line the width is cr = 0.092 ± 0.012 keV. The line 
intensity increases with the 3-10 keV flux, but not as strongly 
as expected in a constant EW scenario, suggesting the pres- 
ence of a constant and a variable Fe Ka line component. 

- The Chandra HETG spectrum has enough energy resolution 
to resolve the narrow component of the Fe Ka line (cr = 
0.027+^-^i8 keV; line intensity (1.5 ±0.2)xlO-5 ph cm'^ s'^). 
The width of the narrow component of the line suggests an 
origin at around 0.2-0.5 pc (~ 30000 rg) from the BH. This 
value is of the same order of magnitude as the molecular 
sublimation radius, suggesting that the narrow component of 
the Fe Ka line might be produced as reflection from the inner 
walls of the molecular torus. If so, because of light travelling 
effects, the intensity of this component has to be constant 
on years time-scales. We assume the presence of a constant 
narrow Fe Ka line (as observed by Chandra HETG) and add 
a second, resolved (now observed to be broader cr = 0.22 ± 
0.04 keV), Fe Ka component with EW= 42+^ eV. There is 
excess emission at 7.06 keV, consistent with being produced 
(at least in part) by the associated Fe KJS emission. 

- For the first time reverberation mapping of the resolved com- 
ponent of the Fe Ka line on timescales of several days-years 
was successfully performed. The intensity of the resolved 
Fe Ka component shows a significant (~ 4cr) 1-to-l cor- 
relation with the 3-10 keV flux variability, however the EW 
stays constant during the 9 years XMM-Newton observed the 
source. The robustness of this result is confirmed by the re- 



sults of the rms spectra which, in a model independent way, 
show an excess of variability at E= 6.45 ± 0.08 keV. This ex- 
cess of variability is consistent with being the resolved com- 
ponent of the Fe Ka line (crS 0.2 keV) varying in such a way 
as to keep a constant EW= 71 ± 36 eV. No measurable lag of 
the reflected component is observed. 

- The width of the resolved component of the Fe Ka line sug- 
gests an origin between 300 and 1000 rg from the BH. This 
location is consistent with the observed Fe Ka variability on 
days to a week timescale and the lack of measurable lag. The 
lack of a relativistic red wing of the Fe Ka line suggests an 
inner radius for the line production larger than several tens 
ofrg (~ 40 rg). 

- The EW= 42^4 eV of the resolved Fe Ka line suggests a 
larger covering factor of the primary X-ray sources (assumed 
to have altitudes of few rg above the BH) compared to the 
one expected from a flat disc annulus, indicating a possible 
azimuthal distribution above the disc of the reflecting mate- 
rial. A possibility is that the material producing the resolved 
Fe Ka emission might be in the form of clouds, perhaps as- 
sociated to the inner BLR (see Costantini et al. 2012). This 
geometry is further reinforced by the consistency between 
the Fe Ka line width (cr = 0.22 keV) and those from the 
broadest components of the UV broad emission lines (Kriss 
et al. 2011). We also observe that the location of the rever- 
berating Fe Ka emission does not correspond to that of any 
X-ray or UV absorption components (Detmers et al. 2011; 
Kaastra et al. 2012; Kriss et al. 2011; 2012; Ebrero et al. 
2011). 

- Significant, but weak (15-20 eV) ionised Fe K emission is 
observed. The ionised emission can be fit equally well with 
two narrow emission lines (from both Fe XXV and Fe XX Vl), 
possibly from a photo-ionised or collisionally ionised gas, or 
by a single broad relativistic emission line (either Fe XXV 
or Fe XXVI). We note that the source high Eddington ratio 
suggests the presence of a standard thin a-disc down to a 
few rg from the BH. However, the neutral Fe Ka line has no 
redshifted wing with no neutral emission closer than ~ 40 
rg from the BH. This suggests that the surface of the in- 
ner accretion disc in Mrk 509 might be highly ionised. For 
these reasons, although the two interpretations for the origin 
of the ionised Fe K emission are equivalent on a statistical 
ground, we slightly prefer the latter interpretation on physi- 
cal grounds. The picture of an higher ionised disc in the inner 
few tens of rg from the BH and less ionised outside is in line 
with the presence of a compact hard X-ray corona, provid- 
ing there a high flux of hard X-ray photons, and a soft more 
extended one, as proposed by Petrucci et al. (2012). 

- A highly ionised, medium outflow velocity (v ~ 0.048 ± 
0.013 c) Fe K absorption component detected in previous 
observations (EW= -13^2 9 ~ 4 cr significance) ap- 
pears much weaker (EW= -3.2^1^ eV, ~ 4 times weaker), if 
not absent, during the 2009 campaign. 

- Previous XMM-Newton observations showed evidence for 
highly ionised high outflow velocity (v ~ 0.05 - 0.2 c) ab- 
sorbers based on a total exposure of ~ 300 ks. We find no 
convincing (> 3cr) evidence for these features during the 
2009 XMM-Newton long (600 ks) monitoring campaign. 
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2009 campaign - summed spectrum - 1 component of the Fe Ka line 
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Table 1. Best fit results of the summed EPIC-pn spectrum of the 10 XMM-Newton observations performed during the 2009 campaign, o: assuming crpe = ctfq Ka',t- assuming 
c^Fe Kfi = o-fq Ka and Npe K/3 = 0.15x Npe K^; f: in units of 10"^ ph cm~^ s~^ 

- Model (1) Single FeKa: Power law + GauspeKa 

- Model (2) Single FeKa + Power law + GauspeK« + GauspeK^e 

- Model (3) Photo-ionised gas + Single FeKa + j3: Power law + GauspeKa + GauspeK^e + Gauspexxv + Gauspexxvi 

- Model (4) Broad line^^xxw + Single FeKa + j3: Power law + GauspeKa + GauspeK^ + Disklinepexxv 

- Model (5) Broad line^exxvi + Single FeKa + /3: Power law + GauspeKor + GauspeK^ + Disklinepexxvi 



2009 campaign - summed spectrum - 2 components of the Fe Ka line 
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Table 2. f: in units of 10"^ ph cm"^ s"^. 



2009 campaign - summed spectrum - 2 components of the Fe Ka line 
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Table 3. f: in units of 10 ^ ph cm ^ s ^ 



- Model (6) Photo-ionised gas + Double FeKa + Absorption line: Power law + GauspeKaN + GauspeKaR + Gauspei^N + Gauspei^R + Gauspexxv + Gauspexxvi - AbsGaus 



- Model (7) Broad ionised Gauss ianFexxwi + Double FeKa: Power law + GauspeK^N + GauspeKaR + GauspeK/5N + GauspeK^R + Gauspexxv + Gauspexxvi - AbsGaus 

- Model (8) lonised-disc-lineYoxxNi + Double FeKa : Power law + GauspeKaN + GauspeKaR + GauspeK^N + GauspejgjR + Gauspexxv + Gauspexxvi - AbsGaus 



